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We present an extension of the corner transfer matrix renormalisation group (CTMRG) method to 
O(n) invariant models, with particular interest in the self-avoiding walk class of models (0{n = 0)). 
The method is illustrated using an interacting self-avoiding walk model. Based on the efficiency 
and versatility when compared to other available numerical methods, we present CTMRG as the 
method of choice for two-dimensional self-avoiding walk problems. 



The self-avoiding walk class of models on the two- 
dimensional square lattice, with a variety of possible in- 
teractions, has mobilised the scientific community for 
about half a century [l|, |2j. The number of exact re- 
sults for such models is limited, and numerical studies 
are hard. A clear illustration of the numerical difficulty 
is the disagreement which existed over the numerical de- 
termination of critical temperature and exponents for the 
standard 0-point model, see for example references 0. 

To date the numerical methods available for the study 
of interacting self-avoiding walks in two dimensions are 
series expansions of walks of lengths of a few tens of 
steps transfer matrices for lattice widths up to about 
12 51 and increasingly complicated Monte-Carlo simula- 
tion methodsQ, limited in practice to only a portion of 
the phase diagram. 

Motivated by these numerical difficulties, we decided 
to extend the Corner Transfer Matrix Renormalisation 
Group (CTMRG) methodQ. The CTMRG method 
is based on White's Density Matrix Renormalisation 
Group method (DMRG)|| and Baxter's corner matrix 
formalism|2| . To date the CTMRG method has only been 
applied to discrete spin models, where it is shown to be 
computationally efficient . 

Our extension to interacting self-avoiding walk models 
exploits the connection betweeen these models and the 
0(n) invariant spin models [l(Tj|. which contains as special 
cases the Ising model (n — 1), the XY model (n = 2) and 
the Heisenberg model (n — 3). The method therefore 
has applications well beyond the self-avoiding walk type 
models (n = 0). 

The O(n) spin model is defined through the partition 
function I 111 
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where (i, j) refers to a sum over nearest neighbour spins. 
The spin Si has n components, and is normed such that 
sf = 1. Another formulation of 0(n) invariant models, 
with the same critical behaviour, is: 
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where the spins are now placed on the lattice bonds [Io|. 
A diagrammatic expansion of Equation [3 follows if we 
identify the 1 as the weight of an empty bond between the 
sites i and j and the K as the weight of an occupied bond. 
This expansion may be expressed in terms of graphs Q 
of non-intersecting loops (collisions at sites are however 
allowed) |10| . The partition function may then be written: 
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where / is the number of loops and b is the number of 
occupied bonds. 

The parameter n is now a fugacity controlling the num- 
ber of loops, and need no longer be taken as an integer. 
This fugacity corresponds to a long ranged interaction, 
since the loops may be of any size. This non-locality is 
undesirable for our purposes. We would like to express 
n as a product over local weights. This may be achieved 
as follows. Each loop may be followed clockwise or anti- 
clockwise. The loops in Equation [3] are not oriented, but 
may be oriented by associating 2^ e ' oriented graphs with 
each non-oriented graph. By associating a loop fugac- 
ity n + (fi-) with the clockwise (anticlockwise) oriented 
loops, the partition function may be rewritten[12(: 
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where Q' is the set of oriented loop graphs, and l + 
(Z_) is the number of clockwise (anticlockwise) oriented 
loops. Setting 71+ = exp(i0) and n_ = exp(— iff) 
gives n = 2cos(#). The oriented loop factor is now 
broken up into local weights by associating a corner 
weight Wi = exp(i#/4) with every clockwise corner and 
Wi = cxp(— i6/4) with every anticlockwise corner. On 
the square lattice, there must be four more corners with 
one orientation, compared with the other orientation, in 
order to close a loop. The product of the local weights 
will then give the correct weight for the oriented loops. 

The partition function may now be rewritten in terms 
of a vertex model llOl: 
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FIG. 1: The 19 allowed vertices in the most general O(n) 
model. K is the step fugacity and r = exp(— /?e), where e is 
the attractive monomer- monomer interaction energy. 



where Vi is the weight of the vertex at site i. The deriva- 
tion given here is only for the simplest case, but we may 
freely change the weights of the vertex configurations in 
order to generate different interactions in the original 
model (see figure 

The limit n — > corresponds to the self-avoiding walk 
modeljl^. In the generalised form presented here this 
corresponds to an interacting self-avoiding walk model 
due to Blote and Nienhuisflfl H3 |. A step fugacity K 
and an attractive short ranged interaction e < are in- 
troduced. In the standard 9 point model the interactions 
are between non-consecutively visited nearest-neighbour 
sites and a given site may only be visited once0,|3|- In 
our current models, this last constraint is relaxed; the 
walk may collide at a site, but not cross, and remains 
self- avoiding for the bonds. The interaction is now as- 
signed to doubly visited sites. An additional weight p is 
added for sites which are visited by a straight section of 
walk (i.e. do not sit on a corner). The partition function 
may then be written: 
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where r = exp(— /3e), Nj is the number of site collisions 
and N c are the number of corners in the walk. This model 
gives rise to the vertex weights shown in figure ^ The 
standard self-avoiding walk model is found setting r = 
and p = 1 and when p = the model has the same critical 
behaviour as the standard 9 point model^^l- During the 
remainder of this letter we shall illustrate the CTMRG 
method for 0(n) models using this Blote-Nienhuis walk 
model. 

Following Baxter0, the partition function of a two di- 
mensional lattice model may be written in terms of the 
product of four matrices representing the four quarters 
of the lattice. The inputs and outputs of the matrices 



are the configurations at the seams of the four quarters. 
These matrices are known as corner transfer matrices. 
In general the four matrices are different, but may of- 
ten be related by lattice symmetries. For our model the 
four matrices are the same up to a complex conjugation 
operation^] . 

It is usually not possible to calculate explicitly these 
matrices for systems with a large number of sites. This 
is where the CTMRG method comes in; the matrices 
for larger lattices are calculated from smaller lattices 
iteratively0- This is done as follows. An initial system, 
consisting of a small number of sites, is mapped exactly 
onto a prototype system made up of four m-state spins. 
At each iteration the system is enlarged by adding sites, 
this enlarged system is then projected back onto the pro- 
totype system in some optimal way, so as to minimise 
the loss of information. The value of m determines the 
amount of information which may be carried forward at 
each iteration, the larger the value of m the better the 
approximation. For details see pRIlfij]. 

As with the DMRG method, the innermost sites of the 
lattice are treated exactly. This means that the CTMRG 
method is most appropriate for the calculations of one- 
point functions such as the site free energy, the site den- 
sity of monomers, specific heats etc. In figure |21 we show 
the density, p, as a function of the step-fugacity K for 
fixed values of r. We chose to fix p = since for this 
value the Nienhuis-Blote model has the same critical be- 
haviour as the standard 9 point model^^l- When r < 2 
the model is expected to be in the self-avoiding walk uni- 
versality class and present a critical transition where the 
density changes continuously from p = for K < K c to 
p > for K > K c . When r = 2 the transition is ex- 
pected to be tricritical, in the same universality class as 
the 9 point, with K c = 1/2 exactly. For r > 2 the tran- 
sition is expected to be first order. These three types of 
behaviour may be clearly seen in figure [2] 

Whilst the raw finite-size data presented in figure [5] al- 
ready gives a fairly precise idea of the behaviour of the 
system (order of the transition and first estimate of K c ), 
practical calculations of critical temperature and critical 
exponent estimates rely on finite-size analysis. In trans- 
fer matrix calculations one usually uses a phenomeno- 
logical renormalisation group based on the correlation 
lengths due to Nightingale^ ■ In the CTMRG method 
an effective transfer matrix may be determined directly. 
However, as already mentioned, the highest precision is 
obtained for the one point correlation functions. Since 
the number of sizes available is an order of magnitude 
larger than that available for transfer matricies it is ad- 
vantageous to exploit the finite size scaling laws for the 
one point functions. In particular it is expected that the 
singular part of the density scales as 



p s (K,L) = L^- 2 p(\K - K C \L X ' V ). 
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FIG. 3: <pl,l-2 as a function of K for the self-avoiding walk 
model (p = 1, t = 0) for L = 10,20,30,40 and 50. The 
horizontal and vertical lines give the corresponding finite size 
estimates of K c and v. 
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FIG. 2: The density, p, as a function of K with p = and 
(A) t = 1, (B) t = 2 and (C) r = 3 for L = 62, 122 and 
182. Estimates of K c from finite size scaling are given in 
table [I] The value of K for the first order transition when 
r = 3 estimated directly from the jump in the density was 
found to be K* = 0.422 ± 0.001 (compared to 0.421 ± 0.001 
found previously 131) • 



This scaling behaviour implies that the function 
\og( Ps (K,L)/ Ps (K,L')) 
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takes the value <Pl,L' = 1/v — 2 when K = K c , indepen- 
dantly of L and L'. Naturally there are additional finite 
size corrections which should be taken into account, but 
the conclusion is that if the function <pl,l'(K) is plotted 
for various values of L and L' then it will converge to a 
fixed point given by f(K c ) = 1/v — 2. In what follows 
we have set L' = L — 2 and looked for solutions of the 
equation 



<p L<L -2( K c) = <PL-2,L-4( K i 
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where K^, is the estimated critical temperature from sizes 
L, L — 2 and L — 4. If such a solution does exist then 
K c = lim^oo and v = lim^oo l/(^ L , L _ 2 (if c L ) - 2). 

Figure 01 shows tp plotted for several values of L for 
the pure self-avoiding walk model (p = 1, t = 0). The 
different curves cross at a point defining K c and v. The 
corresponding estimates of K c and v are shown in Fig- 
ure 0] plotted as a function of 1/L. The extrapolations 
of ICf? and v L to L — > oo are given in Table [I] along with 
prclimcnary estimates for p — and different values of 
t. The accuracy with which the critical points are de- 
termined is in general an order of magnitude better than 
found with transfer matricies. Full details and definitive 
estimates will be given elsewhere 

In (numerically) exact methods, such as the transfer 
matrix method or the exact enumeration method, the 
main source of uncertainty in the results is due to the ex- 
trapolation of a small number of points. In the CTMRG 
method the extrapolation problems are largely removed 
since lattice sizes an order of magnitude larger have been 
reached. The uncertainty now lies in the precision related 
to the calculation of each density point. 
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TABLE I: Estimates of K c from transfer matrix (TM) calcula- 
tions taken from reference [T^jl for p = and from reference^] 
for the self-avoiding walk (p = 1 and r = 0) are compared to 
our (previsional) estimates using the CTMRG. First estimates 
of v using CTMRG are also given. 
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may easily be extended to more complicated interacting 
self-avoiding walk models, such as the hydrogen-bonding 
self-avoiding walk^^ . 
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FIG. 4: (A) Critical point and (B) critical exponent estimates 
for the self-avoiding walk (p = 1, r = 0) as a function of 
1/L. The extrapolated values of K c and v are shown. In (A) 
the three horizontal lines show the previous transfer matrix 
estimate for K c along with the corresponding error bars0. 
In (B) the horizontal dotted line indicates the exact value of 
v = 3/4. 

In this article results for the Blote-Nienhuis interact- 
ing self-avoiding walk are presented since good quality 
results already exist, providing a good test of the effi- 
ciency of the CTMRG method. In particular we have 
focused on the self-avoiding walk model (p = 1, r = 0) 
for which the asymptotic limit is well described by trans- 
fer matrix calculations and yet we still find a substantial 
increase in precision (see Table [|J. It is clear that in 
circumstances where larger system sizes are required to 
extract the scaling behaviour, CTMRG should far exceed 
the numerically exact methods in performance. 

The quality of results presented in this article is virtu- 
ally unattainable with such ease by any other numerical 
method we know of (there is a very small number of ex- 
ceptional cases where better accuracy was obtained 01) 
and so we present CTMRG as the method of choice 
for two-dimensional self-avoiding walk models. CTMRG 
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